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Abstract. Within the framework of the mode-coupling theory of super-cooled 
liquids, we investigate new phenomena in colloidal systems on approach to their glass 
transitions. When the inter-particle potential contains an attractive part, besides the 
usual repulsive hard core, two intersecting liquid-glass transition lines appear, one of 
which extends to low densities, while the other one, at high densities, shows a re-entrant 
behaviour. In the glassy region a new type of transition appears between two different 
types of glasses. The complex phenomenology can be described in terms of higher order 
glass transition singularities. The various glass phases are characterised by means of 
their viscoelastic properties. The glass driven by attractions has been associated to 
particle gels, and the other glass is the well known repulsive colloidal glass. These 
correspondences, in associations with the new predictions of glassy behaviour mean 
that such phenomena may be expected in colloidal systems with, for example, strong 
depletion or other short-ranged attractive potentials. 
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1. Introduction 

The experimental study of tlie glass transition in colloidal systems has been a 
very important test case to assess the validity of theories concerning the formation 
of an amorphous sohd from super-cooled liquids. In particular the mode-coupling 
theory (MCT) has been useful when applied to colloids, modeled as spherical particles 
interacting through a hard sphere purely repulsive potential [H. In this case the MCT 
predicts the existence of a critical volume fraction (f) where the system undergoes an 
ergodic-nonergodic transition, which was observed experimentally using quasi-elastic 
dynamic light scattering j^]. From the physical point of view the MCT describes in a 
fairly accurate way the so-called cage effect, i.e. the fact that at high densities molecular 
motions of particles are constrained by the presence of the surrounding ones which form 
a cage around each particle. Prior to reaching the MCT transition, we may think of 
the particle vibrating within their cages at short time-scales, and escaping from their 
cages at somewhat longer time-scales. The two time scales show up in a well-defined way, 
when the liquid system gets closer to the critical threshold, as distinct relaxation regimes 
separated by a plateau region which becomes more and more extended as criticality is 
approached. Close and above the plateau, there is a power law and the plateau ends 
with another power law prior to enter into the so-called full a regime. The power-law 
region is called f3 correlator in the framework of MCT. The a decay is quite well 
phenomenologically described by a stretched exponential. 

We have briefly sketched the main features of the MCT for colloids considered as 
hard spheres and proceed to consider the novel consequences of adding an attractive 
contribution to the inter-particle potential. In this paper we mainly consider the 
effects of using a hard core followed by a square well potential in order to mimic 
the interactions between colloidal particles. In real systems this can be obtained, for 
example, by covering the surface of the particles with a polymer coating or with depletion 
interaction [^]. A description of such systems using the MCT approach has revealed the 
existence of a set of new and interesting phenomena, that we briefly summarise in what 
follows 0, 0]. In the temperature T and volume fraction plane MCT predicts two lines 
of transition from liquid to glass. One of the lines extends to high temperatures, it can 
be traced essentially to the repulsive part of the potential and tends asymptotically for 
high temperatures to the value (p corresponding to the critical value for hard spheres. 
We will call it the repulsive glass transition line (RGL). On lowering the temperature the 
transition line moves toward higher values of and gives rise to a re-entrant behaviour, 
i.e. one can pass from a liquid to a glassy state either by lowering or by raising T. In 
other words the liquid phase tends to exist in a region which extends more deeply into 
the glassy phase compared to the hard spheres case. The origin of this unusual effect is 
that if the interaction is short ranged enough, there can, in this temperature and density 
regime, be partial cancellation of the repulsive and attractive interactions leading to a 
phenomenon for glasses, not unlike that of a theta point for polymers. 

The other glass transition line, originating in a relatively well-defined energy scale of 
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the well-depth, is almost parallel to the axis and extends on one side of the binodal to 
the other, until it eventually crosses the repulsive transition line. We call it the attractive 
glass transition line (AGL). For sufficiently narrow well- widths and on the high end 
the AGL enters the glass region and terminates in an endpoint. Thus the two sides of 
the line separate two different types of glassy system that we will characterise through 
their mechanical properties. The endpoint of the AGL is a higher order glass transition 
point (called ^3), which will be shown to have special dynamical properties, namely 
the relaxation slows down dramatically on approaching it. Upon increasing the width 
of the potential well, the length of the AGL shrinks and its extension into the glass 
region reduces to zero, and as a consequence the endpoint turns into a higher transition 
point (called A4). We will characterise this rather complex behaviour of the system, 
where repulsion and attraction compete, using two typical dynamical quantities, the 
shear viscosity, which characterises the liquid phase, and the elastic shear modulus for 
the amorphous glass. 

Aspects of the MCT have been tested in many different cases, both experimentally 
and with computer molecular dynamics 0. The most extensive and accurate 
experimental check has been performed in the liquid-glass transition of colloidal systems, 
treated as hard sphere systems, and studied with dynamic light scattering 0. The 
agreement with MCT is quantitatively satisfactory Q]. In a similar fashion computer 
simulations have been used to study the glass transition in simple model systems. 



as diverse as Lennard- Jones binary systems |p, the SPC/E model for water |T0 
orthoterphenil |lT[ and silica [0. In all these cases evidence shows that these systems 
undergo a kinetic glass transition, and the molecular dynamics is well accounted for by 
the idealised MCT of super-cooled liquids. On the other hand there is extensive evidence, 
mainly due to experimental results in colloidal systems, that cannot be interpreted in 
terms only of hard-core potentials. 

Dense systems of colloidal particles characterised by a hard core and strong short- 
ranged attractions have been realized experimentally by adding polymers to either a 
suspension of colloidal hard spheres , in solutions of sterically stabilised particles when 
decreasing the solvent quality |T^, and in copolymer micellar systems when changing 



the temperature Such systems were also studied by Monte Carlo simulations flS 
The facts that cannot be simply explained using hard-core potentials are the following. 

(i) An amorphous material can be formed by increasing the attraction strength 
even though the volume fraction is kept well below the value of the hard sphere glass 
transition [0, [l^, |19|. 



(ii) In mixtures of colloids and polymers, melting of the glass states is obtained by 
increasing the strength of a short range attraction by the addition of small polymers 

(iii) Using solvents of decreasing quality |TB[ in polymer coated colloidal particles 
the long time limit of the density time correlation function at small wave vectors are 
much larger than in hard-sphere systems. 

(iv) Viscoelastic measurements for intermediate frequencies found strongly 



concentration dependent elastic moduli ||I6|, [T^, |18|, |1£ 
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(v) A system in which an anomalous dynamical behaviour has been reported is a 
polymer solution (called L64) where structural arrest accompanied by a logarithmic-like 



decay of the density correlator is observed [20|. 



The paper is organised as follows. In Section ^ we report the details of the 
calculation of the structure factor for a square- well system (SWS), the only input needed 
for the mode-coupling dynamics, which is described in Section ^. The resulting phase 
diagram is described in Section the time dependent density correlation functions in 
Section |^ and finally the various phases are characterised by their mechanical properties 
in Section |. In Section ^ we report our conclusions. 

2. Calculation of the structure factor 

In the framework of the MCT the only quantity that is needed in order to calculate 
the dynamical properties of a super-cooled liquid approaching the glass transition is the 
structure factor. It is well known that, given an inter-particle potential satisfying the 
Ornstein-Zernike equation for the space-dependent pair correlation function g{r), one 
needs a closure approximation in order to solve the equation. Many approximations of 



this type have been proposed and solved in the case of simple model potential In 
the case of the square-well potential the Wiener-Hopf method as formulated by Baxter 
2^ can be used. The interaction potential V{r) for particles a distance r apart, is 



obtained with a hard-core repulsion for r < d, and the negative attractive value — Mq 
within the range d < r < sD. 

The structure factor is specified by three control parameters, i.e. the packing 
fraction related to the hard cores, the temperature T, and the relative width 
e = 1 — d/D oi the attraction shell. In the Baxter approach the structure factor 
Sq is expressed in terms of the Fourier transform Q{q) of the factor function Q{r), 
a continuous real function defined for r > 0, 

=Q{q)Q%q), (1) 

Q{q) = 1 — 2np / dr exp{iqr)Q{r) . (2) 
Jo 



The funcion Q{r) is related to the direct correlation function c(r) |2T| and both 
functionss vainsh beyond the distance D. For < r < D, the following equation for 
Q{r) holds 

rc(r) = -Q'(r) + 27rp ds Q'{s)Q{s - r) . (3) 
In addition for r > 

rR 

rh{r) = -Q'{r) + 27ip ds{r ^ s)h{\r ~ s\)Q{s) . (4) 







For the SWS, g{r) = is fulfilled for < r < d, and therefore, using h{r) = g{r) — 1, 
equation @ splits into three sub-equations. The result for the middle part, A < r < rf. 
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where A = D — d, is very simple since the formula known from the theory for the 
hard-sphere system is reproduced 

Q'{r) = ar + b, (5) 

where the coefficients a and b are introduced by 

rd+A rd+A 

a = 1 — 2-Kp / ds Q{s) , b = 27rp / ds s Q{s) . (6) 
JO Jo 

Defining G{r) = rg{r) one finds, for small distances, < r < A, where A = D — d, 

Q'(r) =ar + b-27rp j'^^^ ds G{s - r)Q{s) , (7) 

Jr+d 

and for the attraction shell, d < r < D, one obtains 



g'(r) =ar + b- G{r) + 2np f ds G{r - s)Q{s) . 

Jo 



(8) 
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Figure 1. The Baxter factor function Q{r) a,t e — 0.03, cj) — 0.504449 and different 
temperatures. The hard-core diameter is chosen as the unit of length. 

A closure approximation for c(r) must be introduced into equation (O) in order 



to complete the system of equations 
approximation |21|, according to which 



and 



c(r) = g{r) 



1 — exp 



V(r)' 



Using the Percus-Yevick (PY) 



(9) 



inside the hard core in equation @ and equation (H) leads to the approximation valid 
ioT d<r < d + A 
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Figure 2. The structure factors function in the PY approximation at various 
temperatures, e = 0.03 and (p = 0.504449. The inset shows the maxima of the structure 
factors. 



G{r) =ar + b-2np f^^^ ds Q\s)Q{s - r) 

J r 

+ 2np r ds G{r - s)Q{s) . (10) 



Equations (H) and (|10D for Q{r) and G{r) are solved numerically choosing on each of the 
three intervals of the variable a grid of equally spaced points Vn, where n = 1, 2, . . . , 1000. 
The integral in equation (H) is then performed to obtain Q{q) and hence Sq. The typical 
values for the factor function Q{r) are shown in figure 0, while figure |^ shows the 
corresponding structure factors Sq for various temperatures, given in term of the well 
depth Mo as khT/uQ. 

In a first attempt to solve the problem associated with the attractive part of the 
inter-particle potential the Baxter potential has been used i.e. the limit of a square- 
well in which the range of the potential vanishes while its depth increases in such a way 
that their product remains constant ||^. The mean spherical approximation has also 
been used in conjunction with an attractive Yukawa potential in an attempt to explain 
colloidal gelation |2^. We may note that more recently other methods of generating the 
structure factors have been explored The main results reported here are found 

not to depend on these details, and we therefore here only comment on the square well 
potential. 
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3. The mode-coupling equations 

The MCT of the ideal liquid-glass transition is capable of describing the ergodic to 
non-ergodic transition in the normalised density correlation function 

In the long time limit 4>q{t) shows a discontinuous variation from a vanishing value in a 
super-cooled liquid to a finite value different from zero fg in a glassy state, on varying 
the external control parameters. In the latter state the system experiences a structural 
arrest. In real systems the transition point is never reached since new dynamical 
mechanisms set in which bypass the ergodicity breakdown. The MCT equations are 

0,(t) + fij0,(t) + uj,{t) + Ql 1^ m,{t - t')4>,{t')dt' = (12) 

where the mode-coupling memory functional ruq is given by 

1 r d^k 
(2^ 

The mode-coupling vertices are determined by the density p, the structure factor Sq and 
the direct correlation function Cg = (1 — l/Sq)/p 

P 



= o / 7^^q,kfkfi^_ki ■ (13) 



% = SgSkS^._,:^- [q-kck + q-{q-k) C|._,^|J . (14) 



The two quantities Qq and Uq are respectively the characteristic frequency of the phonon- 
type motions of the fluid, and a term that describes instantaneous damping arising from 
the fast contribute to the memory function. They are defined as 

mS{q) 
Vq = viq^ 

and z/q = 1 in our calculations. In the long time limit t — > oo, the density correlators 
(pq{t) tend to a value 

. _ (p;(0)P.(oo)) .... 
(IP.(O)P) ^''^ 



the non-ergodicity factor, or Debye- Waller factor. The MCT equations (O) in the static 
limit give rise to the bifurcation relation 

f 1 r d^k 

It is clear that /g = is a solution of equations (ITBp and it corresponds to an ergodic 
state of the system in which the correlations decay for long time. The correlators tend 
instead to a finite value different from zero if the system is kinetically arrested. This loss 



Mode- Coupling Theory of Colloids with Short-range Attractions 



8 




Figure 3. The phase diagram of a square- weU system in the PY approximation for 
a square well relative width e = 0.03. The inset shows the complete phase diagram, 
where the AGL and the RGL are given for various values of e. From j^. 



of ergodicity for is interpreted as the transition to a kinetic glassy state. Therefore 
for some critical values of the thermodynamic parameters, density and temperature in 
our case, bifurcations of the solutions of the asymptotic equations appear that produce 
non-zero solutions. The bifurcations can be multiple, up to the number of control 
parameters of the system. Thus, when a bifurcation gives rise to more than two solutions 
of equations (]TB|), there will exist multiple solutions with finite non-ergodicity factors. 
In these cases one speaks of A^-type bifurcations, with A; = 2, 3, 4, ... In this case, MCT 
predicts that only the state corresponding to the largest value of fq is a stable solution 
of the equations . 

4. Phase diagram 

The phase diagram resulting from the solution of the MCT equations in the large time 
limit is shown in figure ^ We report the liquid-glass transition lines for some values 
of the width parameter e. At this point it is worth reviewing the points made earlier 
about the experimental observations in the introduction. Indeed, the calculated phase 
diagram exhibits the phenomena of points (i) and (ii), and in later sections we shall 
show that all the other points are satisfactorily described. We point out the salient 
features of the phase diagram below. 

(i) For high enough temperatures the repulsive glass lines, RGL, tend to converge 
to a value approaching the critical volume fraction for hard-sphere systems (pc ~ 0.516, 
corresponding to the fact that for large T the short-range part of the inter-particle 
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potential weakly affects the dynamics. In this case the microscopic effect described by 
the MCT is the excluded volume cage effect, the impossibility of the molecules to move 
freely due to the presence of the neighbouring particles. 

(ii) On lowering the temperature a liquid stabilising effect due to the attractive 
forces sets in, so that the liquid-glass lines tend to extend to larger values of the volume 
fraction. This behaviour gives also rise to the possibility of the glass melting when 
lowering the temperature. 

(iii) The reentrant behaviour characterises the interplay of the repulsive and the 
attractive inter-particle forces; the system can become a glass both when lowering or 
raising the temperature. 

(iv) For even lower values of the temperature another phase line, the AGL, appears 
and it is almost parallel to the axis. In this case the particle motions are difficult 
since, due to the attractive interactions, the particle tend to form bonds. 

(v) The two sections of the phase lines, RGL and AGL, intersect at a finite angle for 
small values of the parameter e. The AGL extends for higher volume fractions beyond 
the crossing point and finally ends in a point. 

(vi) The points on the glass lines we have discussed so far all correspond to 
bifurcations of type A2, while the end point of the AGL relates to a higher order 
bifurcation ^3. 

(vii) The extension of the AGL that develops in the glass region leads to the 
coexistence of two different types of glasses. In order to represent the differences 
between the two glasses we present in figure |^ the shape of the non-ergodicity factor fg 
crossing the glass-glass transition. The results are now consistent with the experimental 
observation (iii) in the introduction. The width of the fg is related to the localisation 
length of the particles in the glass. If the width is larger the particles are more localised. 
It is evident that the colloidal particles are more localised in the attractive glass then 
in the repulsive one. It is also possible to show that the localisation of the length in 
the repulsive glass remains more or less unchanged decreasing the temperature, whereas 
for the attractive one the particles become more localised. This seems to confirm the 
idea that the repulsive glass is dominated by the so-called cage effect: when the system 
gets to high packing fraction the particles start to be blocked by their neighbours, 
at certain critical packing fraction 0c each particle is trapped in a cage formed by the 
surrounding particles. Indeed this is a purely geometrical effect and it should not depend 
on temperature. 

(viii) On increasing the relative square well width parameter e the glass-glass line 
tends to shrink and finally coincides with the point of intersection of the RGL and 
AGL when e ~ 0.04. This end-point is a higher bifurcation point of type ^4. We 
note that there is no difference in the structure factor across this curve. Only the late 
stage dynamics, reflected in the values of the non-ergodicity factors, may be used to 
differentiate between these glasses. 
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Figure 4. The typical non-ergodicity factors fq crossing the glass-glass transition line 
for e = 0.03 and cj) = 0.544052. 

5. The intermediate scattering function 

We have studied the dynamical behaviour of the system in different regions of the phase 
diagram, in particular in the proximity of the higher order bifurcation points and 
A4. We will mainly comment on the case where the well width parameter is e = 0.03. 
The wave vectors are given in units of the particles diameter. In figure |^ we present 
different intermediate scattering functions for q = 21.75, obtained by solving MCT 
equations, approaching the glass transition at = 0.504449 and Tc ~ 0.8000. The 
relevant cut across the attractive glass curve is shown in figure |^ using crosses, and 
a transition point. For temperatures well above the glass transition, i.e. T = 0.85, 
4>qit) presents the typical liquid behaviour: after a first short-time decaying due to 
microscopic dynamics, the correlation function starts to relax toward ergodicity and no 
sign of a critical slowing down is evident. When the system becomes close enough to 
the transition temperature (T = 0.80010 in figure ^ the system starts to freeze and 
the typical two step relaxation scenario starts to be evident. We note that we are here 
crossing that transition curve which we believe corresponds to the process of gelation 
at high density. We see, therefore, that we may expect the characteristic behaviour of a 
glass transition in the slow dynamics near gelation at high density. Once the system is 
below the transition temperature (the case T = 0.80000 in figure |), 4>qit) corresponding 
to the gel does not relax to zero anymore, the system is non-ergodic. In this situation 
there is only part of the /3-relaxation toward a non-zero constant value, i.e. the non- 
ergodicity parameter /g 7^ 0. Perhaps it is useful now to discuss some of the MCT 
predictions for the asymptotic behaviour of the intermediate scattering function (see for 
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Figure 5. Time dependence of the density correlators for q = 21.75, e = 0.03, 
(j) = 0.504449 and different temperatures. 

example and references therein). In the time region where t^ '^t <^Ta the so-called 
factorisation theorem holds 

Ut) = n<l) + h,G{t) (17) 

to being a typical macroscopic time, Tq, the a relaxation time scale and the critical non 
ergodicity parameter. Equation |1^ implies a universal behaviour in which wave- vectors 
and time factorize. G{t) is a scaling function which describes the whole relaxation 
pattern. Near the glass transition the behaviour of the function G{t) may be calculated 
analytically in the proximity of the a and (3 relaxation. It can be shown that the 
asymptotic behaviour can be expressed in terms of a rescaled time r = t/tc, where to- 
is a time scale that depends crucially on the distance to the transition, i.e. on the 
separation parameter a 0. Close to the glass transition, for T > T^, G{t) is given by 



G{t < 1) ~ (18) 
G{t > 1) ~ -T^ (19) 

where the two exponents are related to the exponent parameter A, which is obtained 
from a stability matrix [|], via the relation \ = T'^{l-a)/T{l-2a) =T^{l + h)/T{l + 2h), 
which implies < a < 1/2 and < 6 < 1. The scaling behaviour of the correlation 
function changes in the proximity of higher bifurcation points. In particular it is possible 
to show that near an ^3 point G{t) ~ l/ln(t/ti) whereas near an it goes like 



G{t) ~ l/ln^(t/ti), where ti is a microscopic time scale . This is quite a remarkable 



behaviour and represents a rather strong prediction of the theory. There have been 
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reports of such behaviour in some system and it will be interesting to see if this 

phenomenon is discovered in a variety of different colloidal systems. 



6. Mechanical properties 

Two particularly interesting quantities to investigate in the vicinity of the glass 
transition are the shear viscosity and the elastic shear modulus. Those two quantities 
can be measured in real systems and such measurements can provide a good insight in 
the glass transition phenomena and its relations with MCT. The complex shear viscosity 
r]*{(f), T, uj) can be calculated in terms of the normalised time correlation function of the 
density fluctuations (pkit) [H, Q 

1 2 



60 n'^Jo 



die 



iuit 



dkk" 



and is related to the complex shear modulus G 

G*{(j),T,uj) = i UJ ri*{(j),T,u). 
It is possible to take the ut ^ limit in equation pO 



dlnSk 
dk 

T, uj) by 



60 vr^ 



dkk' 



dlnSk 
dk 



1 2 



fk 



(20) 



(21) 



(22) 



and this gives the static shear viscosity. We note that this formula, in that it depends 
on the square of the correlation function, is different in structure from that used by 
Weitz and co-workers, where the stress-relaxation function is linearly related to the 



density correlation function ||2^. This difference would lead, in principle, to substantially 
different predictions for the loss modulus, and in particular around its minimum. It 
seems likely that the idea the stress modulus is linear in the correlation function is 
more suited to visco-elastic properties dominated by large domains, and their changing 
surface area under shear |2^, whilst the dependence on the square of the correlation 
function is more suited to dense fluids where the properties are dominated by loss of 
mobility due to caging effects. All the quantities in the equations (pOf), ( pTD and (p^) 
can be evaluated for the present system, and therefore it is possible to use the numerical 
solution of equation (^) in order to solve them. This has been accomplished performing 
the integrals with standard numerical integration over 2000 g-vector varying between 
and 500. In figure ^ the static shear modulus has been represented for three different 
packing fractions {(f) = 0.539672, = 0.544052 and = 0.580000) as a function of 
the temperature. These three cases correspond to crossing the glass-glass line, crossing 
the A3 point and a density where the repulsive and attractive glasses have become 
indistinguishable. In the first case it is clearly possible to distinguish between the two 
glasses. For low temperatures there is strong dependence of the elastic viscosity on the 
temperature, in particular the system becomes more and more rigid on decreasing the 
temperature. When the system crosses the glass-glass transition there is a discontinuity 
in the elastic response which clearly indicates that the structure is changed. Increasing 
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Figure 6. The static shear moduH on crossing the glass-glass transition line for 
e = 0.03 and close to the A3 bifurcation point. 



further the temperature the elastic behaviour does not change so much anymore and 
for large temperature the system behaves like a hard spheres suspension. In this case 
the glass is originated by the cage effect and consequently the particle are forced to 
move inside a fixed volume that does not change with temperature. It is then evident 
how the differences between the two glasses are manifested by mechanical properties. 
We have noted that the difference in the shear modulus as the point is approached 
is described a by a power law, and it is also possible to show that in this regime the 
differences in mechanical properties of the two systems are due to the different shape 
of the non-ergodicity parameter fg, and hence the long time residual motions in the 
gel, and not to the contribution of the equilibrium structure factor |]29|. We have 



performed this calculation also for the dynamical quantities defined by (|20D and (|21|). It 
is evident that, if the system is in the glassy phase we have for the long time behaviour 
of the density correlators (f)g(t 00) = /g 7^ 0, and consequently the integral in (^) 
is unbounded. Thus the solution posses a zero frequency pole. In order to obtain 
convergent solutions for a glass we have replaced (pq(t) with 4>q(t) = (pg(t) — fg which 
decays to zero at infinite time. In figure |^ and figure 0b the real and the imaginary 
part of the complex viscosity, rj'{w) and ri"{u), are presented for a fixed volume fraction, 
(f) = 0.504449, for different temperatures, approaching and crossing the glass transition. 
In figures ||a and |b G'^u) and G"{u) are shown for the same points as figure |^. Since 
figures |^ and ^ contain the same information we discuss only figure |[ G'{uj) represents 
the elastic response of the material: the higher its value the stiffer the material at 
that frequency scale. The imaginary part of the shear modulus describes the viscous 
behaviour of the system and so the dissipation. Both quantities can be easily measured 
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Figure 7. Relaxation of the shear viscosity on crossing the attractive transition Hne 
for e = 0.03 and cj) = 0.504449. 



experimentally, for example as the response of a sample to small oscillatory shears which 
weakly perturb the system ||2^. Approaching the glass transition we note that the G'{lS) 
begins to develop a shoulder at low frequency. This effect is due to the slowing down of 
the dynamics of the system, in other words the formation of the plateau in the is 
responsible for the formation of a region where G'iuS) varies only slightly with frequency. 
The range of such a region tends to increase approaching Tg . Indeed a similar behaviour 
has been observed experimentally in measurements of linear viscoelasticity in a colloidal 
suspension with hard sphere interaction ||2^ . If the system is in a glassy state the value 
of the real part of the static shear modulus G'{u) = 0) is finite, indicating that the 
system is solid and consequently it presents an elastic behaviour. In the liquid phase, 
however, the system does not show any elastic behaviour so the G'{uj = 0) tends to zero. 
It is then clear that at the glass transition the system abruptly changes its behaviour, 
presenting a singularity in the static shear modulus. In figure ^ the behaviour of the 
imaginary part of the shear modulus is represented. At high frequency the curves for 
different temperature are the same, showing the microscopic dynamics which is the 
same for all the temperatures. Approaching the glass transition a second maximum 
starts to emerge. Such a maximum represents the a-relaxation and it moves toward 
low frequency on decreasing the distance to the glass transition temperature Tg. The 
minimum of G"{uj) corresponds to the plateau region in time, and the power law behavior 
in frequency on both sides has been already observed in hard spheres systems p7 |. 
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Figure 8. Relaxation of the shear modulus on crossing the attractive transition line 
for e = 0.03 and ^ = 0.504449. 



7. Conclusions 

In this paper we have outhned some of the experimental observations that have 
been made of particle gels in colloidal systems where we know there to be a strong 
short-ranged potential. We have shown that all of these features may be interpreted 
within the paradigm that such particle gels represent a new type of attractive glass. 
However, having made this correspondence, we then find a number of predictions of 
new phenomena that have yet to be confirmed experimentally, or are just being so. 
These include the re-entrant behaviour for the phase diagram, the logarithmic dynamic 
behaviour near this re-entrant regime, and in the general the characteristic alpha 
relaxational behaviour near gelation. We have calculated the mechanical manifestations 
of these phenomena also. Thus we have shown how the glass-glass transition may be 
studied using the shear modulus, and that the particle gel must be expected to be 
much stiffer than the colloidal glass, but that the differences between them vanishes 
near the A3 point. We have also studied, to our knowledge for the first time, the 
frequency dependent modulus and loss modulus for these systems, using the paradigm 
that the system is undergoing a glass transition. MCT type treatments seem to be 
suited to study colloidal systems of this type because, at least of the repulsive system, 
the glass transition seems to lack some of the fluctuations of molecular hquids. We also 
can calculate phase diagrams, and mechanical and scattering properties from a single 
theory. We expect that, in future, dense particle gels will be better fltted to the theory 
present here, than previous phenomenological treatments. 

There remains much to be done in many directions. It now begins to seems likely 
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that the glass paradigm is suitable for interpretation of particle gelation. However we 
must await further detailed experimental work to see just how extensive the agreement 
with MCT type treatments will be. We may also speculate that this type of idea will 
have much broader significance that particle gels. There is every reason to hope that 
polymer gelation and protein gelation may also find a description in similar terms. 
Certainly this is a rich area for future exploration. 
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